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Abstract 

Quantum field model of unstable particles with random mass is suggested to de- 
scribe the finite-width effects in decay rate. Within the framework of this model we 
derive the convolution formula for a width of the channels with unstable particle in 

in ■ 

a final state. The distribution function of random mass is considered for unstable 

<n : 

■ particles of arbitrary type. 
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1. Introduction 

Quantum field description of the unstable particles (UP) with a large width runs into some 
problems, which are under considerable discussions pp. These problems have both the con- 
ceptual and technological status and arise due to UP lie somewhat outside the traditional 
formulation of quantum field theory [2] . Unstable field can not be treated as asymptotic state 
and perturbative approach is unfit in the resonance neighborhood. This conceptual prob- 
lems is connected with methodological difficulties, such as ambiguity in definition of mass 
and width. Therefore, the new quantum field approach [2] (Bohm et al), phenomenological 
models [3] and effective theories of UP jlj are actual now. 

Convolution method [5] is convenient and clear phenomenological way to evaluate the 
instability or finite-width effects (FWE). This method describes FWE in the processes of 
type $ — > 0i0 — > 010203-.., where is UP with a large width. The intermediate unstable 
state is simulated by the final state in the decay $ — > 0i0 with invariant mass, described 
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by Breit-Wigner-like (Lorentzian) distribution function. The phenomenological expression 
for a decay rate has convolution form 



where p(q) = Mr^(g)/7r|P(g)| 2 . In Eq.Q p(q) is probability density of invariant mass 
distribution, P(q) = q 2 — M 2 + iMV^q) (Altarelli et al), T($ 010(g)) and T^q) 
are partial width of $ and total width of in the stable particle approximation, when 
mi = q 2 . The formula for a decay rate, which has a close analogy to the Eq. ([T]) . was applied 
first to describe FWE in B and A decay channels with p(770) and ai(1260) in the final 
states, which have large total widths [SJ. It was shown that the contribution of FWE to 
the decay rates of these channels are large (20-30 %) and that the account of it significantly 
improves a conformity of experimental data and theoretical predictions. Analogous results 
were obtained in Ref. |3] for the dominant decay channels of $(1020), p(770) and A'* (892). 
The decay rates of the near-threshold decay channels t — > WZb, cWW, cZZ were calculated 
with help of convolution formula (CF) in Ref. 0. It was shown in these works, that the 
FWE play a significant role in the near-threshold processes. 

The convolution formula was derived by direct calculation from the decay-chain 
method in Ref. [Zj, where in analogy with inclusive processes the contribution of all de- 
cay channels of UP is described by function p(q 2 ) = qF(q) /\P(q)\ 2 ''. The essential elements 
of this derivation for vector and spinor UP was the expressions r] mn = —g mn + qmln/l 2 and 
f] = q + q as numerators of vector and spinor propagators (q = q^)- The convolution 
formula was derived for the decay chain t —>■ bW — > bfifj in the limit of massless fermions 
/ [S](Galderon and Lopez-Castro). Quantitative analysis of convolution and decay-chain 
calculations of the t — > WZb decay rate was fulfilled in Ref. [3] (Altarelli et al). The formula 
for a decay rate, which is in close analogy with was received in Ref. [2] for the scalar 
UP within the framework of the "random mass" model. An UP is described in this model 
as quantum field with a "smeared" (fuzzy) random mass in accordance to uncertainty prin- 
ciple for energy and lifetime of unstable quantum system 8J. The FWE is connected with 
fundamental principle, which gives the relation 5m *r « 1, that is 5m w V in the rest frame 
of reference (5E = 5m, c = h = 1) 3 J. So, uncertainty principle leads to the interpretation 
of kinematic value q 2 in Eq.(^) as random mass square. Thus, the intermediate states of 
UP, which are traditionally defined as virtual, in the neighborhood of q 2 = M 2 are not differ 
from real ones in accordance to uncertainty principle. This interpretation is connected with 
a smearing of mass shell and with above mentioned definition of r] mn and r), which are pro- 
portional to the polarization matrix for vector and spinor UP (see section 3). As was noted 




(1) 
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in [7] , this proportionality leads to the factorization of expression for width in decay-chain 
method, and, as consequence, to the CF ((Tj). Thus, the suggested model is theoretical basis 
of the convolution method on the ground of uncertainty principle. 

In this paper we consider the generalization of the model which includes vector 
and spinor fields. Within the framework of this generalized model the CF is derived for 
UP of arbitrary type. To determine the probability density p(m), which is an analogue 
of p(q) in Eq. ([T|) . we put a connection between the model and effective theory of UP with 
modified propagators, used in Ref . (7j . It was shown that this connection leads to Lorentzian 
probability density p{m), which was commonly used in convolution method. Suggested model 
is applicable to the decay processes of type $ — > 4>i4>(q) and gives the convolution formula 

for UP of arbitrary type. 

2. The model of unstable particles with a random mass 

The effect of mass smearing is described by the wave packet with some weight function w(p), 
where /i is random mass parameter H . The model field function, which simulates UP in the 
initial, final or intermediate states, is represented by the expression: 



In Eq.(J2I) $ a (a;,/i) are the components of field function, which are determined in the usual 
way when m? = p is fixed (stable particle approximation). The limits of integration will be 
defined in the sections 3 and 4. 

The model Lagrangian, which determines "free" unstable field has the convolution 

form: 



In Eq.(jnj) /i)) is standard Lagrangian, which describes model "free" field /i) with 

fixed mass m 2 = p. 

From Eq.© and prescription <9$(x, /i)/<9$(x, p) = 5(p — //) it follows Klein-Gordon 
equation for each spectral component: 



As a result we have standard momentum representation of field function for fixed mass 
parameter /i: 




(2) 




(3) 



(□ -n)$ a (x,fj) = 0. 



(4) 
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All standard definitions, relations and frequency expansion take place for <E> a (fc, //), but the 
relation = \Jk 2 + p defines smeared (fuzzy) mass-shell due to random p. 

The expressions (|2J) and (jHJ) define the model "free" unstable field, which really is some 
effective field. This field is formed by interaction of "bare" UP with it's decay channels and 
includes nonperturbative self-energy contribution in the resonant region. Such an interaction 
leads to the spreading (smearing) of mass from p st (p) = 5(p — M 2 ) for the bare particles 
to some smooth density function p(p) = \u(p)\ 2 with mean value p ph M 2 and ph T. So, 
the UP is characterized in the discussed model by the weight function oo(p) or probability 
density p(p) with parameters M and T (or real and imaginary parts of pole). A similar 
approach has been discussed by Matthews and Salam in Ref. jH]- 

The commutative relations for a model operators have an additional ^-function: 

lK(k,p), $+(g,/i')]± = 5{p - p)5{k - q)S aP , (6) 

where subscripts ± correspond to the fermion and boson fields. The presence of 5(p — p) 
in Eq.© means an assumption - the acts of creations and annihilations of particles with 
various p (random mass square) don't interfere. So, the parameter p has the status of 
physically distinguishable value as random m 2 . This assumption directly follows from the 
interpretation of q 2 in Eq. ((!} as random parameter p. By integrating both side of Eq.© 
with weights lu*(p)uj(ij!) one can get standard commutative relations 

fc{k),*tm± = 8(k-q)Safi, (7) 
where is full operator field function in momentum representation: 

= / *i(M)w(A*)dA*- (8) 

It should be noted that Eq.((ZJ) follows from Eq.© when f \u(p)\ 2 dp = 1. 

The expressions (0) and © are the principal elements of the discussed model. The 
weight function uj{p) in Eq.© (or p(/i)) is full characteristic of UP and the relations (0) 
define the structure of the model amplitude and of the transition probability (section 3). 
The probability density p(p) will be defined in the fourth section by matching the model 
propagator to renormalized one. 

With help of traditional method one can get from Eqs.(j2j), (pjj and © the expression for 
the unstable scalar Green function [Sj: 

(0\T(<f>(x),<Ky))\0) = D(x-y)= [ D{x - y, p)p{p)dp . (9) 
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In Eq.(jnj) D(x, /i) is denned in the standard way for the scalar field with m 2 = /i and describes 
UP in an intermediate state: 



-ikx 



D(x,u) = - — - / — dk. (10) 

v ' w (2vr) 4 J k 2 -fi + ie v ' 

The right side of the Eq.fjUJ) is Lehmann-like spectral (on /i) representation of the scalar Green 
function, which describes the propagation of scalar UP. Taking into account the connection 
between scalar and vector Green functions, we can get the Green function of the vector 
unstable field 

Analogously Green function of the spinor unstable field: 

D{X, fl) = (id + y/flD{z, fl) = ^31 / k 2^ u e ~ lkXdk > ( 12 ) 

where k = k^ 1 . These Green functions in momentum representation have a convolution 
structure: 

Dmn(k) = J D mn (k,fj,)p(fj)dfj,, D(k) = J D(k } fi)p(fi)dfj,. (13) 

3. Model amplitude and the convolution formula for a 
decay rate 

In this section we consider the model amplitude for the simplest processes with UP in a final 
state and get the CF (JTJ) as direct consequence of the model. The expression for a scalar 
operator field j^j: 

^0*) = 7^ I "Wv I ^-%^-e^dq, (14) 



(2n)V 2 J vw P J ^Tq 

where g° = \Jq 2 + p and a ± (q, /i) are creation or annihilation operators of UP with momen- 
tum q and mass square m 2 = /i. Taking into account Eq.tJHJ) we can get: 

[a-fcvM+WU [r(x),a + (k tf i)].= e **», h» = y/& + ^. (15) 

(27r) 3/ V 2A: ° 

The expressions ()15|) differ from standard ones by the factor uj(fi) only. From this result it 
follows that, if d + (k, /f)|0) and (0\d~(k,fi) define UP with a mass m 2 = [i and a momentum 
k in the initial and final states, then the amplitude for the decay of type <3> — > (jxf>\ has the 
form: 

A(k,fj) =u(fi)A st (k,n) , (16) 



where A st (k, /i) is amplitude in a stable particle approximation when m 2 = p. This amplitude 
is calculated in a standard way and can include high corrections. Moreover, it can be effective 
amplitude for the processes with hadron participation [HIE!- 

To define the transition probability of the process $ — * 4>4>%, where <fi is UP with a large 
width, we should take into account the status of parameter p as physically distinguishable 
value, which follows from Eq.©. Thus, the amplitude at different p don't interfere and we 
have the convolution structure of differential (on k) probability: 

dT(k) = J dT st (k,ii)\uj(fi)\ 2 dii. (17) 

In Eq. (117)1 the differential probability dT st (k, p) is defined in the standard way (stable particle 
approximation) : 

dT st (k, p) = —S(kg, - JL - k^A^ik, p)\ 2 dksdh , (18) 

where k = (/c$, k^, k\) denotes the momenta of particles. From Eqs. (fT7|) and (fTS|) it directly 
follows the known convolution formula for a decay rate 

r(m$,m 1 )= / T st (m^,mi, p)p(p)dp , (19) 
J Ho 

where p(p) = |o;(/i)| 2 and p , \i m are defined in Refs. E] as threshold and maximal 
invariant mass square of unstable 0. 

An account of high corrections to the amplitude ()16|) and, hence, to Eq. (|19J) keeps convo- 
lution form ()19j) . This form can be destroyed by the interaction between the products of UP 
{(f)) decay and initial $ or final <fti states. The calculation in this case can be fulfilled in a 
standard way, but UP in the intermediate state is described by the model propagator. How- 
ever, a calculation within the framework of perturbative theory (PT) can not be applicable 
to the UP with large width, that is to the short-living particle. In any case, the applicability 
of PT, model approach or convolution method to the discussed decays should be justified 
by experiment. The correspondence of CM to the experimental data was demonstrated for 
some processes E] , but this problem needs in more detailed investigation. In this 
connection we should note the analysis of higher-order corrections for processes with UP 
[1]. The separation between factorizable and non-factorizable corrections make it possible 
to build the effective theory of UP jl] . 

When there are two UP with large widths in a final state $ — > 0102, then in analogy 
with the previous case one can get double convolution formula: 

T(m q> ) = J J r st (m$;n 1 ,[i2)p 1 (n 1 )p 2 (H2)dnidii2- (20) 
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The derivation of CF for the cases when there is vector or spinor UP in a final state can be 
done in analogy with the case of scalar UP. However, in Egs- lfTljl . (jl5j) and (JTfij) we have a 
polarization vector e m (q) or spinor u% (q), where q is on fuzzy mass-shell. As a result we 
get polarization matrix with m? = //. For the vector UP in a final state: 

5^ e ™(<?) e n(<?) = -9mn + 9m9n/ > • (21) 

e 

For the spinor UP in a final state: 

E^'^rfo) = i(?T^- (22) 

In Eqs. fl21|) and (|2*2*|) sum run over polarization and g° = a/<? 2 + [i- 

The formulae ()19|) and (|20|) describe FWE in full analogy with the phenomenological 
convolution method e 5j and with some cases of the decay-chain method [SI El- Thus, we con- 
sider the quantum field basis for CM, which takes into account the fundamental uncertainty 
principle and is in good agreement with experimental date on some decays. To evaluate 
FWE for the case, when UP is in an initial state, we must account the process of UP gen- 
eration. When UP is in an intermediate state, then the description of FWE is equivalent to 
the traditional one, but the model propagators are determined by Eqs.Q - (fT3|). 



4. Determination of p{p) from renormalized propagator 

The possibility of p(/i)-determination directly follows from the connection of the decay-chain 
method (DCM) and convolution method 7\. As was shawn in Ref. [7j, this connection 
leads to the convolution formula where in accordance with uncertainty principle q 2 is 
interpreted as smeared mass square parameter /i, which distribution is described by the 
expression: 

1 y^r(^) 

In Eq.(J33J) r(//) is /i-dependent full width and P(/i) _1 is propagator's denomenator. It 
should be noted, that the convolution structure of Eq. (0) and universal structure of Eq (j2*H|) 
don't depend on the definition of P(fi). It has a complex pole structure \i — \ir and can 
be approximated by the Breit-Wigner /i — M 2 + iMT([i) j^j or another phenomenological 
approximation. The expression ()23)1 is very simple and convenient in practical calculations 
of decay rate, where the error of approximation is small. 

Here we'll consider the definition of from the matching model propagators to stan- 
dard dressed ones jS] . This consideration is rather methodological than practical and demon- 
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strates the connection between model and traditional descriptions. Let us associate the model 
propagator of scalar unstable field (J£5J) with standard one: 

k 2 - n + it * k 2 - ml - U(k 2 ) ' 1 ' 

where U(k 2 ) is conventional self-energy of scalar field. With help of an analytical continuation 
of the expressions ()24j) on complex plane k 2 — > k 2 ± it and prescription jS] : 

U(k 2 ± it) = ReU(k 2 ) t UmU(k 2 ) (25) 

the conformity ({23$ can be represented by the equality 

- d V= 13 o,J ± ., mn »,n (26) 



o k 2 — jj±it k 2 — m 2 (k 2 ) ± ilmll(k 2 ) 
where m 2 {k 2 ) = m\ + _ReII(£; 2 ). With account of round pole rules and dp = <i(p =F p(p T 
ze) = p(p) =F 0(it) two Eqs.([26p can be combine into the equality (p ± it — > z): 

, dz = I \ (27) 



z — k 2 k 2 — m 2 {k 2 ) — ilmli{k 2 ) k 2 — m 2 {k 2 ) + ilmll{k 2 ) 
The left side of Eq.([27|) is Cauchy integral, which equal to 2mp{k 2 ) and after a change k 2 — > p 
in the final expression for p we have: 

/ N 1 M W f? M 

PlW tt [p - m 2 (p)] 2 + [M(/i)] 2 ' 1 J 

The expression ()28p for p(fc 2 ) in Breit-Wigner approximation is usually exploited within the 
framework of convolution method. From Eq. (j28|) and definition p(p) = |u;(p)| 2 it follows: 

p — m 2 (p) ± ilmll(fi) ^ ^ 



The ambiguity of sign in ()29)1 is not essential because the expression |u;(p)| 2 only enters into 
the physical values. In the parametrization ItoII(p) = v //ir(p) we have relativistic Breit- 
Wigner u;(p) and Lorentzian p(p), which coincides with the expression (|2*3)) for renormalized 
P(q 2 ). Inserting the expression (}2*Hj) into the left side of Eq. (j24j) one can check with help of 
Cauchy method the self-consistency of Eqs. ()24j) and ()28|) . 

Thus, we have put the correspondence between the model |2j - [0] and some effective 
theory of UP with renormalized propagator of scalar UP. To establish such a correspondence 
for the vector UP we insert p(p) into the model propagator ([T5|) with D mn (k,fi), defined by 
(JTTJ) for vector unstable field: 

f 00 -g mn + kmkn/fi 1 Imll(fi) _ 

J k 2 -/i + it tt [p-m 2 (p)] 2 + [Jmll(p)] 2 P ~ 1 j 

~9mn ~\~ krnk-n/ p ^ I 1 ]fip 



2iir J k 2 — jj, + it p — m 2 (p) — ilmll(jj) ji — m 2 (p) + iimll(p) 
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With help of Eq.(|2*5|) and above used method we can represent the second part of Eq. pUj) in 
the form(/i — > z — p, ± ie): 

1 f dz Qmn k m k n /z Qmn k m k n /k (31) 

2m J z — k 2 z — m?(z) — ilmll(z) k 2 — m?{k 2 ) — HmU.{k 2 ) 

The right side of Eq. (pH|) coincides with the expression for propagator of vector UP, which 
leads to the convolution formula (JTJ in the decay-chain method jjj. The numerator of this 
effective propagator coincides with r) mn (k), which was used in [7j. In Eqs.([3*0|) and (|3*T|) the 
value Il(fc 2 ) is defined for vector field as transverse part of polarization matrix PQ. The 
calculations of U(k 2 ) in effective theory (unstable hadrons) or in gauge theory (Z,W-bosons) 
can run into some difficulties. In the first case loop calculation can be ambiguous and we 
should use traditional Breit-Wigner approximation m 2 (/i) m M 2 and ImU(p) « pT(p). To 
escape the gauge- dependence in the second case we can use pole definitions of mass and 
width [T]. 

The description of p(p) by the universal function (|28|) for scalar and vector fields can be 
justified by the general structure of parametrization for bosons: 

m 2 (q 2 ) = m 2 + ReU(q 2 ), ImU(q 2 ) = qY(q 2 ) . (32) 

In the case of unstable fermion we have another parametrization scheme: 

m(q 2 ) = m + ReY>(q 2 ), JmS(g 2 ) = Y(q 2 ) . (33) 

So, we need in additional analysis to define fermion function p(/i). If we choose for fermion UP 
the universal density function (J23j) . which follows from convolution method [7], then we must 
do exchange ImU(fi) —->■ y/JlImYl(fj,) in the Eq.([28|). Inserting the result into Eq. (jl3j) with 
D(x,n), defined by Eq.([T2|). we can get the correspondence between the model propagator 
of fermion unstable field and the effective theory one: 

k + ^JI k + k 

p(ji)dfi — > — , (34) 



k 2 — \x + ie k 2 — m 2 (k 2 ) — zfcS(/c 2 ) 



where k = a/ (kk). The numerator of the right side of Eq. ()34j) coincides with the expression 

The transitions ({211) , an d establish the correspondence between the discussed 
model and some effective theory of UP within the framework of traditional QFT approach. 
These transitions follow from the determination of p(fi), that is from the accounting of 
interaction, which forms the wave packet (0) and mass-smearing. The most important feature 
of the effective theory, chosen in such a way, is the possibility to connect the decay-chain 
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method and convolution method within the framework of this theory jjj. So, we have 
some self-consistency of the discussed model, effective theory, convolution and decay-chain 
method. However, due to some difficulties, which arise in traditional approach, the search 
of alternative p(/i) - definition is actual now. 

5. Conclusion 

The finite width effects in the processes with participation of UP can be described by renor- 
malized propagator, decay-chain method, convolution method and effective theory of UP. 
The convolution formula is convenient instrument for calculations of decay rate and gives 
the results in accordance with experiment. In this paper we have considered the model of 
UP with a random mass and derived the convolution formula as a direct consequence of the 
model. The model operator function and Lagrangian have a convolution structure, which 
describes mass-smearing in accordance with uncertainty principle. 

The principal element of suggested model is probability density function p(p), which 
describes the main properties of UP. Traditional description of UP in the intermediate state 
by resonance line with complex pole (or by dressed propagator with mass and width as 
parameters) corresponds to the model description of UP in arbitrary state by function p(p) 
with the same parameters. We have considered the determination of p(p) from DCM and 
by matching the model propagator to renormalized one. This approach is equivalent to the 
convolution method or truncated decay-chain method. 

The second p(/i) - determination has some restrictions, caused by propagator renormaliza- 
tion peculiarities. The question arises, also, concern the possibility to describe mass-smearing 
of bosons and fermions by the universal function p(p). Moreover, as the mass-smearing effect 
follows from the fundamental uncertainty principle, then the search of p(/i) from the first 
principles is reasonable. It should be noted also, that the model erases a difference between 
the real and virtual states of UP at peak region. 
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